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We investigate effects of doped holes' hopping on Neel order in the two-dimensional Hubbard 
model. Semiclassical staggered moments are computed by solving saddle point equations derived 
from a path-integral formalism. Effects of quantum fluctuations are taken into account by the 
Schwinger boson mean field theory. We argue that hopping of doped holes is ineffective in suppressing 
Neel order compared to rapid supprestion of Neel order in high-temperature superconductors. After 
destruction of Neel order, the quantum disordered phase sets in. Taking the strong coupling limit in 
the quantum disordered phase leads to a model of spinless fermions and bosons but no gauge field 
interaction. 
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I. INTRODUCTION 



In high-temperature superconductors, one remarkable 
feature is rapid destruction of Neel order by hole doping. 
In fact, only 2% doping hole concentration is enought to 
suppress Neel order. While critical disorder is 50% for 
the bond percolation threshold and 41% for the site per- 
colation threshold. How do doped holes suppress Neel 
order in such an effective way? 

Naive expectation is that the hopping process of doped 
holes suppresses Neel oder. In this paper, we examine the 
effect of disorder brought by hopping of doped holes. As a 
model, we take the Hubbard model because at half-filling 
it is reduced to the S = 1/2 antiferromagnetic Heisen- 
berg model that describes the undoped parent compound 
of high-temperature superconductors, and perhaps it is 
the simplest model to see disorder effect by doped holes' 
hopping on Neel order. We use path-integral formal- 
ism developed by Schulz. [1,2] We solve saddle point 
equations to compute the magnitude of the staggered 
moment. Since saddle point equations are semiclassical 
equations, solutions do not contation effects of quatum 
fluctuations. We take into account quantum fluctuation 
effects in terms of Schwinger boson mean field theory. [3] 

The remainder of this paper is organized as follows: In 
sec. II, we rewrite the Hubbard model following Schulz. 
[1,2] In sec. Ill, we dervie saddle point equations. Solv- 
ing these equations, we compute doping dependence of 
the semiclassical staggered moment. We examine effects 
of quantum fluctuations on Neel order by the Schwinger 
boson mean field theory. We show that there is the quan- 
tum disordered regime in which quantum fluctuations 
suppress Neel order. In sec. IV, we discuss the effective 
action in the quantum disordered phase. Sec. V is de- 
voted to summary and discussion. 



II. PATH INTEGRAL FORMULATION OF THE 
HUBBARD MODEL 

In order to investigate doped holes' hopping effects on 
Neel order, we first rewrite the model in a convenient 
form following Schulz. [1,2] In the coherent state path- 
integral formulation, the partition function of the model 
is given by Z = JVcDcexp{-S) with S = S + Su, 
where 



S = 



13 



(It 



'J 



Sn = 



3 



(1) 



(2) 



Hereafter r dependence of fields is implicit. The 
summation * s t & ken over the nearest neigh- 

bor sites. Carrier fields are represented in a 
spinor: c t = T ( Cj| ) and Cj = ( Ci\ Cii ) ■ 
Using the identity, [4] the on-site Coulomb inter- 
action term can be rewritten as, U^2i-n^riji = 

-(U/4) Ej [Kt + n 3l? + {cj<rcj) 2 + (n iT + n i; )] , 
where the components of the vector er = (a x , a y , a z ) 
are the Pauli spin matrices. Introducing Hubbard- 
Stratonovich fields for the charge and spin fluctuations, 
[1,2] we obtain Z = J VcVcVSVcj) c cxp(-So - Su), 
where. 





Su — I dr 
lo 



) c j Cj Cj 



(3) 
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up to constant. Here the vector Sj represents the lo- 
calized spin moment. The scalar <p c j is associated with 
charge fluctuations. For the charge degrees of free- 
dom, we take the uniform value at the saddle point: 
<t> C j = (cjcr^Cj) = 1 — 5, with 5 the doped hole concentra- 
tion. 

Thus, the approximate action is given by 



S = S 



dr 



CjfTCj + 



U^S 



(4) 



Note that the first term in the square brackets has the 
form of Hund coupling between the localized spin mo- 
ment and the carrier's spin. 



III. EVOLUTION OF NEEL ORDER BY HOLE 
DOPING 

Now we consider doped holes' hopping effects on Neel 
order. Our strategy is the following: First, we esti- 
mate semiclassical staggered magnetic moments by solv- 
ing saddle point equations. Secondly, we examine sta- 
bility of Neel order against quantum fluctuations. For 
analysis of quantum fluctuations, we use Schwinger bo- 
son mean field theory. [3] 

In the Neel ordered phase, the localized spin moment 
Sj has the following form: 



s, = (-iyse z . 



(5) 



Substituting this into Eq.4, and then performing Fourier 
transformations, we obtain 



S = ^2 ^2 ( c k (iuj n ) c k+Q (iuj n ) ) 



-ioj n + efe - n -USc 
—USa z —iu r , 
+/3NUS 2 , 



C k+Q(iu n ) 



(6) 



where = — 2i(cos 

kx cos ky ) and the summation in k- 
space is taken over half of the first Brillouin-Zone. The 
energy dispersion of the carriers is given by ±£k with 



Ek = y e 2 . + {US) 2 . After integrating out fermions, we 



obtain 



-/3(E k -n) 



* = - 2 E'{ ln [i 

k 

+ ln[l + e-^- £fc -^]}+/3A^5' 2 (7) 

At zero temperature, variation of the action with re- 
spect to the chemical potential \x yields 



1 ]T [6{E k + n) + 9(-E k + M )] = 1 - S 



(8) 



and variation with respect to S yields 

iEj [9{Ek + ^ ~ e{ ~ Ek + M)1 = h (9) 
k k 

where 9{x) = 1 for x > and zero otherwise. 

We compute S by solving saddle point equations (8) 
and (9). Figure 1 shows U/t dependence of S at half- 
filling. The same result was obtained in Ref. [5] based on 
the Hartree-Fock approximation. Figure 2 shows doping 
dependence of S at U/t = 6, 8, 12. Figure 3 shows effec- 
tive exchange interaction J(x)/J = (2S*) 2 . Temperature 
dependence of the antifcrromagnctic correlation length 
is given by £af — 0.26 exp(1.38J(a;)/T) according to the 
renormalization group analysis of non- linear sigma model 
and numerical simulations. [6] 

Note that non-zero values of 5* do not necessarily im- 
ply that there is Neel order. Because solutions of saddle 
point equations neglect quantum fluctuation effect. We 
need to examine stability of Neel order against quantum 
fluctuations. For this purpose, we apply the Schwinger 
boson mean field theory. [3] In the Schwinger boson mean 
field theory, the localized spin moment Sj is represented 

by 



Sj — 2 Z ' jCrz ^ 



(10) 



with the constratint X) ct =t l z 3 aZ i a ~ We- derive the 
effective action of the localized spin system by integrat- 
ing out the carrier fields. Detail of calculations is given 
in Appendix A. The result is, 



-"spin 



/ 

JO 



dr 



^2 Zj a d T Zj a AijAij 



(id) 



(11) 



Zi±Zji and 



where J = 4t 2 /U and Aij = z^Zji 
z^Zji — Zi±Zjf. (In terms of the fields Sj, S sp i n is, 

•Sspin = J2j Sf my + J* drJ Y,(i,j) S» • Sj, where Sf my 
denotes the Berry phase term for the localized spin mo- 
ment Sj. This is nothing but the spin S antiferromag- 
nctic Hcisenberg model.) 

In the Schwinger boson theory, Neel order is stabilized 
if S is larger than S c = 0.19660 as shown in Ref. [3]. 
(This result is briefly summarized in Appendix B.) The 
parameter regions of S > S c and < S < S c are shown 
in Fig. 4. In the S > S c regime, Neel order is stabilized. 
Whereas in the < S < S c regime, Neel order is sup- 
pressed by quantum fluctuations. The point S = S c can 
be taken as a quantum critical point. [7] The < S < S c 
regime is called the quantum disordered regime. [8] Note 
that this quantum disordered regime is identical to that 
in Ref. [8]. The condition of S > S c for the stability of 
Neel order turns out to be g < g c in Ref. [8]. 

In Fig. 4, the 5 = regime is also shown. Contrary to 
the other regimes, there is no antifcrromagnctic Heisen- 
berg type correlations in this regime. Therefore, spin 
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fluctuations in this regime are disconnected to the origi- 
nal spin correlations at half-filling. 



Note that the hopping for the fcrmions from j-site to 
i-site containes a matrix, 



IV. QUANTUM DISORDERED PHASE IN THE 
STRONG COUPLING LIMIT 



UiUi 



F ■ 



_Aij 
F- 



(16) 



In this section, we discuss the effective action for the 
quantum disordered regime in the strong coupling limit. 
The effective action in the quantum disordered regime is 
given by 



Jo 



dr 



3 (id) 



u 



3 

+ Yz ja d T z ja +J2 X i : 



3° 



J 



3° 



(id) 



(12) 



This action consists of the free fermion part and the 
Schwinger boson part that describes the spin S antifer- 
romagnctic Hcisenberg model. The interaction between 
them is of the form of Hund coupling. The action (12) is 
so-called the spin-fermion model. 

Now let us consider the strong coupling limit, U/t — > 
00. Since the Hund coupling term is dominant in this 
regime, we move to a frame in which the Hund coupling 
term is diagonalized. In such a frame, the doped hole's 
spin at the j-site is in the direction of the localized spin 
moment Sj. The transformation to this frame is given 

by 

(13) 
(14) 



c j ~ Ujfj, 

where the matrix Uj is given by 



After this transformation, the action reads 



dr 



-tYUiUiUifj + fjUjUifi) 

{id) 
3 

^ Zj<yd T Zj a + (^ZjcZja _ ^) 

3° ja 

- 2 Y (zjizif - z 3 \Zi\) {zifZji - znzjt) 

(id) 



(15) 



with Fij = z^Zj^ +ZiiZji and A i3 = z^z^ - z^Zj-f. 
Note that F^ describes ferromagnetic correlations and 
Aij describes antiferromagnetic correlations in the local- 
ized spin system. In fact, (Fij) is taken for the mean 
field in the ferromagnetic spin system and (A^ ) is taken 
for the mean field in the antiferromagnetic spin system 
in the Schwinger boson mean field theory. [3]. 

A similar action can be derived in the slave-fermion 
mean field theory of the t-J model. However, there is a 
crucial difference. If we take the strong coupling limit of 
U/t — > 00, then one finds that the fermion hopping only 
couple to ferromagnetic correlations in the localized spin 
system. Therefore, coupling between fermions and the 
gauge field that describes antiferromagnetic fluctuations 
is absent. 



V. SUMMARY AND DISCUSSION 



In this paper, we investigate effects of doped holes' 
hopping on Neel order. What we have found is that dis- 
order effects induced by holes' hopping on Neel order is 
rather small. In fact, the critical doping concentration is 
S c ~ 0.40 at U/t = 10. This value is substantially larger 
than that in high-temperature superconductors. There- 
fore, hopping processes of doped holes are not so effective 
in suppressing Neel order. For destruction of Neel order 
in high-temperature superconductors, holes must behave 
like an excitation which suppresses Neel order more effec- 
tively. Such an excitation would be intimately connected 
with properties of the localized spin system. 

After destruction of Neel order, the quantum disor- 
dered phase appears. (This quantum disordered phase 
is special to the antiferromagnetic correlations. In fact, 
there is no such phase in case of the ferromagnetic cor- 
relations because there is no quantum fluctuations that 
suppress the semiclassical long-rage order as in the anti- 
ferromagnetic case.) If we take the strong coupling limit 
U/t — * 00, then the coulping between the doped holes and 
the antiferromagnetic fluctuations are lost. Here we first 
derive the spin-fermion model starting from the Hubbard 
model. After representing the localized spin moments 
by the Schwinger bosons, we take the strong coupling 
limit. There is another way of taking this strong cou- 
pling limit. If we take the strong coupling limit first at 
the spin-fermion model, it is believed that the model is 
reduced to the t-J model. Applying the slave-boson the- 
ory, we obtain a system of spinless fermions and bosons 
with a gauge field interaction. In the derivation of the 
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slave-boson representation of the t-J model, one big as- 
sumption is that there is a deconfincment phase of a U(l) 
gauge field theory. Whether there is a dcconfinement 
phase or not is still an unsolved issue. By contrast, there 
is no need to assume a deconfincment phase when we 
introduce boson fields to describe the localized spin mo- 
ments at the spin-fermion model. After taking the strong 
coupling limit, we obtain a system of spinless fermions 
and bosons but there is no gauge field interaction. Our 
analysis suggests that taking the strong coupling limit is 
not justified or indirectly suggests that there is no decon- 
fincment phase. A situation such that taking the strong 
coupling limit is not allowed occurs when there is a term 
of spin-orbit coupling like H so = iY^a j) c\\ij ■ crcj +h.c. 
In the presence of such a term, doped holes rotate their 
spin at every hopping process. [9] Therefore, we cannot 
take the strong coupling limit. 



We integrate out / and ff 



'off 



-Trln 



d T — /i + Ujd T U. 



U 



(A4) 
,• = 



where tij — t for the nearest neighbor sites and U 
otherwise. 

We expand the logarithm in Eq. (A4) with respect to 
tn. The second order term is 



0(2) _ i T 



1 



u tij Ui Uj 



d T - ii + Uid T Ui 
1 

d r — [I + Ujd T Uj — H^O, 



~tjiUj Ui 



(A5) 



Applying the derivative expansion technique, we obtain 
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APPENDIX A: DERIVATION OF THE 
ANTIFERROMAGNETIC HEISENBERG MODEL 

In this appendix, we derive the antifcrromagnetic 
Heisenberg model from So + Su, where So is defined by 
Eq. (1) and Su is defined by Eq. (4), by applying second 
order perturbation theory with respect to t. [10]. In order 
to describe the localized spin moments Sj, we introduce 
the Schwinger bosons through Eq.10. with the constraint 
S(t=t i z ja z ja = We rotate the carrier's spin in the 
direction of the localized spin moment at the same site 
in terms of the following unitary transformation: 



where 



z n 



z j1 



(A2) 



The total action reads 



S= [ dr 
Jo 



^fjidr-^ + UjdrU^fj 



-t^ifiUiUjfi+fjUjUifi] 

(id) 



1 



1 



iuJ n (ij) (7,(7 

■0 



* lUJ n + fl+ |(T iuj n + fl + ^(J 1 



X / dT{<T\U i U j {T)\o'){o'\U i Ui{T)\o) 

Jo 

+(higher derivatives) (A6) 

After the summation over the fcrmion Matsubara fre- 
quencies, we take (3U — > oo limit. Thus, we obtain 

S^ = ~ fdr^AiiAj, (A7) 

<*>j> 

where J = At 2 /U and — z^Zj± — z^z^ and Aij = 

ZilZjl -ZilZjl- 

On the other hand, the expansion of the logarithm in 
Eq. (A4) with respect to Ujd T Uj gives the Berry phase 
term for the localized spin moments as follows. The term 
with the first order of Ujd T Uj is, 



S T" = -Tr- — 



UjdrUj. 



(A8) 



(Al) Applying the derivative expansion technique, we obtain 



P tZ a n + M + "J CT j Jo 

(A9) 

After the summation over the fermion Matsubara fre- 
quency, we take (5U — > oo limit. Thus, we obtain 



rrBcrry 
^eff 



J ° 3° 



(A10) 



This is nothing but the Schwinger boson representation 
of the spin's Berry phase. In terms of the original spin 



(A3) moment fields Sj, the action S< 



-tBcrry 



spin — «-> e ff 



out to be the antifcrromagnetic Heisenberg model. 



(2) 

S cS turns 
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APPENDIX B: THE COMPUTATION OF S c 

Here we briefly summarize the result of Ref. [3] for 
S c = 0.19660. 

In terms of the Schwinger boson fields, which is defined 
in Eq. (10), the antiferromagnetic Heisenberg Hamilto- 
nian H = J^(i j) &i ■ Sj, reads 



We introduce mean fields A-,- = 



t J 



( Z it Z 



A 3 - A* - A 

system is given by 



(Bl) 

z i\ z n - z a z 3\) and 

Zjf) and assume the uniform value 
const. Then, the free energy of the 



+ J A 2 - A(2S + 1), 



(B2) 



where A is a Lagrange multiplier to impose the con- 
straint J^cr z jcr z jv = 2S 1 , and w k = ^/ \ 2 — A A 2 J 2 a 2 , with 
ak = (sin k x + sin k y )/2. 

The variation with respect to A and A yields 



IV^coth^ * 



k 



1 1 
— > — coth 



(B3) 



(3uj k 2S + 1 



k 



2 A 

At zero temperature, Eq.(B4) has the following form: 



(B4) 
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FIG. 1. The localized spin moment S versus U/t at 
half-filling. 



with K the complete elliptic function and p — 2AJ/X. 
The right hand side is a monotonically increasing func- 
tion with respect to p and it takes the maximum value 
S c = 0.19660 at p = 1. Therefore, Eq. (B5) has a solution 
of p < 1 for S < S c . Solutions of p < 1 imply that spin 
wave excitation have gap, and there is no Bose-Einstein 
condensation of Schwinger bosons. [11]. Therefore, there 
is no Neel order for S < S c . 
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FIG. 2. The semiclassical staggered moment S versus the 
doping concentration S for U/t = 6, 8, 12. 
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FIG. 3. The doping dependence of the effective exchange 
interaction J(6)/J = (2S) 2 for U/t = 6, 8, 12. 
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FIG. 4. The Neel order regime S > S c (Neel) and the quan- 
tum disordered regime < S < S c (QD) on the U /t-S plane. 
In the S = regime, there is no antiferromagnetic correlation 
that is associated with the Heisenberg antiferromagnetic type 
correlation. 
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